Mode-locking in driven vortex lattices with transverse ac-drive and random pinning 
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We find mode-locking steps in simulated current-voltage characteristics of driven vortex lattices with 
random pinning when an applied ac-current is perpendicular to the dc-current. For low frequencies 
there is mode-locking only above a non-zero threshold ac force amplitude, while for large frequencies 
there is mode-locking for any small ac force. This is consistent with the nature of transverse temporal 
order in the different regimes in the absence of an applied ac-drive. For large frequencies the 
magnitude of the fundamental mode-locked step depends linearly with the ac force amplitude. 
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PACS numbers: 74.60. Ge, 74.40.+k, 05.70. Fh 

Non-linear dynamics of vortices driven by a cur- 
rent in random media leads to several interesting non- 
equilibrium phases, such aSpplastic flow, moving smec- 
tic and moving Bragg glass. Er~l3 These dynamicaLphases 
can be characterized bvj-their temporal orderotjlD u and 
mode-locking responsesJiHT 2 ] When a vortex array with 
average intervortex spacing a is moving at a high enough 
velocity v , it is possible to have temporal order at the 
washboard frequency luq = 2irv/a, which .results in a 
peak at ujq in the voltage power spectcwm.Ercl This has 
been observed in numerical simulationsoll3 and also re- 
cently in experiments £m When the system is driven by 
a dc + ac force with frequency fl, interference phenom- 
ena leads to mode-kicking steps for vortex velocities such 
that ujq = (p/q)fLli2rtL3 This interesting effect has been 
observed experimentally by Fiorytj and by Harris et al.13 
Recently, we have numerically studied how the existence 
of mode-locking in driven vortex lattices depends on thp 
presence of temporal order in each dynamical regime.!! 2 ] 

Mode-locking phenomena has been extensively studied 
in other systems-in the past, e.g., Josephson junctions 
(Shapiro steps) ,£3 Josephson junctipa arrays,EJ super- 
conductors with p eriod ic pinningJUtJ and charge den- 
sity waves (CDW) .113123 Driven vortex lattices with ran- 
dom pinning have two important features that distin- 
guish them from these systems, (i) There is no in- 
herent periodicity, as for example in Josephson junc- 
tion arrays and superconductors with periodic pinning. 
Temporal order and periodicity are induced dynamically 
due to the vortex-vortex interaction, which tends to fa- 
vor a structure, close to a triangular vortex lattice at 
large velocities] 2 ! (ii) The vortex displacements are two- 
dimensional vectors. This is an important difference with 
respect to-,CDW systems where the displacement field is 
a scalar.E2l In particular, the behavior of the displace- 
ments in the direction perpendicular to the driving fp 
shows phenomena like a transverse, .critical current 
and a transverse freezing transitional] at high velocities. 
It can therefore be interesting to study the possibility 
of mode-locking when an ac force is applied in the direc- 
tion perpendicular to the direction of the dc driving force. 
Recently, it has been found in rectangular periodic pin- 



ning arraysEll and in Josephson junction arraysEl that a 
transverse ac force leads to a new type of "transverse" 
phase-locking in these cases. In this paper we will inves- 
tigate the possibility of transverse mode-locking in driven 
vortices with random pinning. __ 
The dynamics of a vortex in position rv is given byO'El 



dr ■ 



>V) 



F(t), (1) 



where = |r; — r 3 -| is the distance between vortices i, j, 
Ti v = |r s ; — r p | is the distance between the vortex i and a 
pinning site at r p , r\ = $ "^ c2d is the Bardeen-Stephen 

friction and F(t) = ^-[Jdc + J ac cos(r2i)] x z is the 
driving force due to an alternating current J ac cos(Sli) 
superimposed to a constant current Jdc- The vortex- 
vortex interaction is considered logarithmic: U v (r) = 
-A„ln(r/A), with A v = $ 2 /8ttA, and A = 2\ 2 /d is 
thtveffective penetration depth of a thin film of thickness 
dBO The vortices interact with a random distribution 
of attractive pinning centers with U p (r) = —A p e~( r /^ , 
£ being the coherence length. Length is normalized by 
£, energy by A v , and time by t = rj^ 2 /A v . We consider 
N v vortices and N p pinning centers in a rectangular box 
of size L x x L y , and the normalized vortex density is 
n v = N v £, 2 j L x L y = Bt; 2 /&o. Moving vortices induce a 
total electric field E = — v x z, with v = Vj. 

We study the response of the vortex lattice to a dc force 
plus a transverse ac force, F = Fd c y + F ac cos(fM)x solv- 
ing Eq. (1) for different values of F ac and f2. The simula- 
tions are at T = for a vortex density n v = 0.04 in a box 
with L x /L y = V3/2, and N v = 64, 100, 144, 196, 256, 400 
(we show results for N v — 256), and we consider weak 
pinning strength of A p /A v = 0.05 with a density of 
pinning centers being n p = 0.08. We impose periodic 
boundary conditions and the long-range interaction is 
determined by Ref. ^Tj. The time integration step is 
At = O.OOlr and averages are evaluated during 131072 
steps after 3000 stepp-for equilibration. 

In a previous worklia we studied the case of a longitudi- 
nal ac force, relating the mode-locking response with the 
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presence of temporal order for the longitudinal compo- 
nent of the velocity. Here we analyze now the transverse 
temporal order from the transverse power voltage spectra 
(corresponding to the transverse velocity), which are 



is mode-locking. We therefore choose the values of Cl 
such that the expected step, V s t ep — Cla/2ir, would cor- 
respond to velocities V belonging to a given dynamical 
regime of the limit F ac — 0. Each simulation is started 
at (v y ) « 0.97557a/27r with an ordered triangular 
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FIG. 1. Velocity- force curve around the main interference 
condition V = Clao/2n for three typical drive frequencies Cl. 
Each case show results for two values of amplitude F ac (the 
curves are shifted in Fdc for clarity). Insets show correspond- 
ing voltage power spectrum for F ac = and V ~ V s t ep . Verti- 
cal dashed line in the spectral density indicates the expected 
washboard frequency u>o- (a) Cl = 0.02, F ac — 0.01 (left), 
F ac = 0.03 (right), (b) Cl = 0.04, F ac = 0.02 (left), F ac = 0.08 
(right), (c) n = 0.19, Fac = 0.09 (left), F ac = 0.23 (right). 

shown in the insets of Fig. 1. They are calculated as 

rT „ Tr ,. s , s , 2 
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FIG. 2. (a) Velocity-force curve around the main interfer- 
ence step for Cl = 0.08 and F ac = 0.3. (b)-(c)-(d) Typ- 
ical time averaged coarse-grained density of vortices for a 
mode-unlocked state bellow the step, a mode-locked state 
in the step, and a mode-unlocked state above the step, re- 
spectively, (c) Typical voltage power spectrum for the three 
ac-driven regimes mentioned above. 

lattice up to values such that (v y ) « 1.02551a/27r by 
slowly increasing the dc force Fd c with AFd c = 0.00005 — 
0.00025. For low Cl, for which we have plastic flow when 
Fac — > 0, we find that there are no interference steps in a 



(F c < F dc < F p , F c 



0.01, F p w 0.03). In this case 



S x (uj) — J Q dtV x (t) exp(iuit)\ 2 at the different dy- 
namical regimes for F ac — 0B The first regime above 
the critical depinning force F c is the plastic flow regime wide range of F ac (shown in Fig. 1(a) for F ac /V step < 1 

(left curve) and F ac /V s t ep > 1 (right curve)). For inter- 
mediate CI, for which we have smectic flow when F ac — > 0, 
we find that there are no steps for small amplitudes, 
FaclV s tep < 1, while there are steps for F ac /V step > 1, 
as shown in Fig. 1(b) in the left and right curves, respec- 
tively. For high CI, corresponding to a transverse solid 
regime when F ac — > we find that there are steps both 
for small F ac /V s t e p < 1 and large F ac /V s t e p > 1 values 
of the ac amplitude, as we can observe in Fig. 1(c). We 
therefore find-a behavior similar to the case of longitudi- 
nal ac forces.113 In the present case, when the dynamical 
regime has transverse temporal order, any small amount 
of F ac will induce transverse mode-locking, while for the 
dynamical regimes that do not have transverse temporal 
order, a non-zero (threshold) value of F ac is needed to 



we find a broad band spectrum without temporal or- 
der [inset of Fig. 1(a)]. Similar behavior is found in 
the "smectic flow" regime (F p < F dc < Ft, F t w 0.06), 
shown in the inset of Fig. 1(b). This is reasonable, since 
we know, that the transverse motion is diffusive in both 
regimes.Q Only for Fd c > F t , in the "transverse solid" 
regime, we find clear evidence of temporal order in the 
transverse velocity. This is seen in the inset of Fig. 1 (c) 
where well developed peaks appear at the washboard fre- 
quency, loq, and its harmonics. We are now ready to 
study the response to a superimposed transverse ac-force 
F ac cos(Clt), for varying values of F ac . For a given Cl, 
we expect the main interference step (p = q = 1) to oc- 
cur when V — V s t ep = Cla /2tt (i.e., Cl = loq) if there 
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induce transverse mode-locking. 

In Fig. 2(a) we show in detail a typical V — Fd c curve 
around the transverse mode-locking step. To visualize 
the spatial structure of trajectories in the transition we 
define a coarse-grained vortex density p v (r,t). We take 
a coarse-graining scale Ar < a,Q. In Fig.2(b)-(c)-(d) we 
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FIG. 3. (a) Low frequency voltage noise in the transverse 
and longitudinal direction around the main interference step. 
The dashed lines indicate the mode-locking transitions, (b) 
Phase of the washboard frequency component of the longitu- 
dinal voltage Fourier transform around the step. 

show the temporal average (p v (r,t)) of the density for 
three typical values of Fd c> corresponding to voltages V < 
V step (Fig. 2(b)), V = V step (Fig.2(c)), and V > V step 
(Fig. 2(d)). We observe in Fig. 2(b) that within the mode- 
locked state vortices follow one-dimensional trajectories. 
The wavy nature of the trajectories is, of course, due to 
the transverse ac force. Figures 2(e)-(f)-(g) show typ- 
ical transverse voltage spectral densities S x (uj) for the 
three cases mentioned above. We see that there is a sig- 
nificant reduction in the width of the washboard p^ftk, 
within the step; a typical signature of mode-locking£3'E2l 
In Fig. 3 (a) we show the low frequency voltage noise 
in both directions, perpendicular P x and longitudinal 
P y to the dc-force, defined as P X)JI = lirn^o S XiV (oj). 
We see that also the low frequency noise is greatly re- 
duced within the step. (There is a noise peak inside the 
step which corresponds to a transition between differ- 
ent mode- locked structures). In Fig. 3(b) we show the 
phase <p(u)o) of the washboard frequency component of 
the longitudinal voltage Fourier transform V(ojq), defined 
as V(ojq) = y/ S y (u>o) exp(i0(ojo)). Here we see explicitly 
that within the "phase-locked" state there is a well de- 
fined "phase" which varies within the range < <f> < tt. 

In Fig. 4(a-b) we show the range (width) AFd c for the 
case F p < F dc < F t and F t < F dc respectively. The error 
bars and the mean values were estimated by repeating 
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FIG. 4. Step width AF dc vs F ac 2w/Qao = F ac /V s t ep - (a) 
Q = 0.04. (b) fi = 0.13 (A) points and fi = 0.19 (O) points. 
Solid line shows a fit to A\ Ji{F ac /V s tepv3)\. 

the simulation for three different disorder realizations. In 
Fig. 4(a) we show AFdc f° r H = 0.04 vs F ac , which cor- 
responds to the smectic flow regime for F ac — ► 0. We 
see that there is mode-locking only above a non-zero 
threshold value F ac /V s tep ~ 1 (see inset of Fig. 1(a)). In 
Fig. 4(b) we show AF dc for two frequencies £1 = 0.13, 0.19 
vs F ac , which correspond to the transverse solid in the 
F ac = limit. We can collapse (approximately) both 
curves into a single curve if we plot AFdc v s F ac /V s tep- 
Our results follow closely a dependence of the form 
AFdc ~ A\Ji(F a c/V s tepV3)\ with A being a constant. 
In the inset of Fig. 4(b) we can see that there is a lin- 
ear dependence of the mode-locking intensity with F ac - 
This is very different frota, [transverse mode-locking in 
periodic pinning systems Ji3't3 in which the step width 
follows AF dc oc {F ac ) 2 . The rather surprising result that 
in the random pinning case the transverse mode-locking 
intensity has a linear i^c-dependence can be explained 
as a consequence of the existence of transverse temporal 
order in the F ac — limit. We can show this with a very 
simple effective model. The moving lattice can be de- 
scribed approximately by an equation of motion for the 
velocity v of its center of mass, 

v = F dc + F ac cos(fif) - ^2 GU g sin ( G ' r ) ■ (2) 

G 

The Uc are the components of an effective periodic force, 
due to the interaction of the nearly periodic moving lat- 
tice (with reciprocal vectors G) with disorder. For weak 
disorder (small Ug) a first order correction can be ob- 
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tained assuming that in zero order r = ro + (v)i + 
F ac sm(Qt)/n. This gives for the instantaneous veloc- 
ity v and average velocity (v) the following expressions 
at first order in F ac : 



v = F dc - G Ug 

G 

G.F. 



sm 



./o 



+ 2Jj 



G.^r + (v)t 
G.F ac 



(v) = F dc - J2 GU g sin(G.ro) | J 

-2j 1 f^V(G-<v>-n)l 



sin(Qt) 



G.F a 



(3) 



<5(G.(v)) 



(4) 



We consider now an anisotropic triangular lattice with 
one of its principal axes parallel to Fd c = yFd c , and 
for simplicity we keep only the shortest reciprocal vec- 
tors {G} = {g s , gi, gi - g s } where g s 



gi = x^-^= + y^. 
£7 gl _g s + AUi and U s , 



x—A? and 
Then, we consider J7 gl = Ui = 
= U s . Here AUi represents a 
small deformation of the perfect triangular lattice in the 
y direction. U s and Ui could be related, respectively, to 
the smectic and longitudinal structure factor peaks of the 
moving vortex system. For F ac = F ac x we obtain for the 
transverse velocity in the limit F ac = 0, 



2itAUl . 
v T = =r- sm 



a 



V3 



271" / I.. 

— (ro + (v)t) 



(5) 



With this approach we obtain the phase-locking range 
for the first interference step = u>o, 



a 



a 



(6) 

Even when it was derived for small F ac and small dis- 
order, Eq. (6) is the approximate relation found in the 
simulation, shown in Fig. 4(b). From equations (5) and 
(6) we see that, for small F ac , temporal order in the trans- 
verse direction is directly related through AUi with the 
linear dependence of AF dc on F ac . It is interesting to 
note that in the perfectly periodic case, AUi = 0, there 
is no "transverse temporal order" and the mode-locked 
step widths would be quadratic in F ac . This is because 
in the periodic case vortices would move in straight lines 
without any transverse component of the velocity, and 
transverse mode-locking would arise as a second order 
effect. In conclusion, a small amount of disorder (ran- 
dom pinning) is enough to induce transverse temporal 
order, and, thus, a linear step width dependence with 
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